Experimental Hamiltonian identification for controlled two-level systems 
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We present a strategy to empirically determine the internal and control Hamiltonians for an 
unknown two-level system (black box) subject to various (piecewise constant) control fields when 
direct readout by measurement is limited to a single, fixed observable. 
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Accurate determination of the dynamics of physical 
systems and their response to external (control) fields 
is crucial for many applications, especially quantum in- 
formation processing (QIP). Quantum process tomogra- 
phy (QPT) is a general procedure to identify the unitary 
(or completely positive) processes acting on a system by 
experimentally determining the expectation values of a 
complete set of observables 

However, QPT is only a partial solution in many cases. 
One problem is the assumption that one can experimen- 
tally determine the expectation values of a complete set 
of observables, or at least perform arbitrary single qubit 
measurements for a register of n qubits. However, most 
QIP proposals rely on a single readout process, i.e. mea- 
surement in one basis. For example, qubits encoded in 
internal electronic states of trapped ions or neutral 
atoms [sj are read out by quantum jump detection via 
a cycling transition 4]. In solid-state systems based on 
Cooper-pair boxes j^, Josephson junctions j^, or elec- 
trons in double- well potentials j7j, final readout is via 
charge localization using single electron transistors @ or 
similar devices. Finally, solid state architectures based 
on electron or nuclear spin qubits @ , are expected to be 
limited to Uz measurements via spin-charge transfer. 

Generally, one assumes that a single projective mea- 
surement per qubit is sufhcient since arbitrary one- 
qubit measurements can be realized by performing a lo- 
cal unitary transformation before measuring to achieve 
a change of basis. However, implementing these ba- 
sis changes requires accurate knowledge of the dynam- 
ics of each individual qubit and its response to control 
fields in the first place. Yet, in particular for solid-state 
qubits d, Ifl S US J it is difficult to predict the actual 
dynamics of a qubit precisely based on computer models 
and theory alone, since they may be sensitive to fabrica- 
tion variance and even vary from one qubit to the next 
as a result. We address this problem of system identifi- 
cation with a single fixed measurement basis and present 
a general strategy to identify the internal and control 
Hamiltonians for a qubit subject a number of piecewise 



constant controls using a single, fixed readout process. 
We can then bootstrap this process to implement QPT. 

Although the technique may be adapted to dissipative 
systems, in this paper we assume Hamiltonian evolution 
for the purpose of identifying the dynamics relevant for 
the implementation of local unitary operations, i.e. de- 
coherence times much greater than the gate operation 
times, which is crucial for systems of interest in QIP. 
The system evolution is thus governed by a unitary op- 
erator U{t,to), via p{t) = U{t,to)p{tQ)U{t,tQy , where 
p{t) is the system density operator, and U{t, to) satisfies 
the Schrodinger equation 



in^U{t,to) ^ H[i{t)]U{t,to), 
at 



(1) 



where H is the Hamiltonian of the system. The case of 
interest is when H depends on external fields {/m}, which 
we can control. Assuming the fields are sufficiently weak 
and act independently, H has the form 
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where Hq is the free evolution Hamiltonian and for 
TO > describes the interaction of the system with field 
/,„. Each H„i can in turn be expanded in terms of the 
Lie algebra generators; for a two-level system, the Pauli 
matrices aj for j € {x,y,z}: 



(3) 



Thus we need to determine the real constants dmx, dmy 
and dmz for to = 0, 1, . . . , M (the dmo can be ignored 
since they result only in an unobservable global phase 
factor) . 

Geometrically, we can represent the state of the sys- 
tem by its Bloch vector s(t) = {sxt Sy, Sz)^ where Sj = 
Tr[(Tjp(t)]. If the system satisfies Eas. (|ll3() . its Bloch 
vector evolves as 
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S(i) = i?o+ V /rn(<)i?m s(t) 
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FIG. 1: Rotation about axis d and observable z{at). The 
expectation value of the measurement operator will oscillate 
sinusoidally, amplitude and frequency are directly related to 
the declination and length of the vector representing the ap- 
plied Hamiltonian. 



where are the real anti-symmetric rotation generators 
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and R = Rq + fmRm generates a rotation about 

(6) 
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d = do + ^ /md„ 



with d„i = {dmx,dmv,d„izV for m = 0,1, . . . , M . If 
fm{t) vary in time then so does d. If they are piecewise 
constant, having fixed values /,„ for ^ ^ ^ii then 

s{t) for t G [ioj^i] rotates about the fixed axis d, where 
|d|| is the rotation frequency, and the unit vector d — 
■jjg^d specifies the rotation axis. This allows us to give an 
explicit formula for the trajectory of s{t) with s(0) = Sq: 

s(at) — So costtf + d(so • d)(l — cosQf) + (sq X d) sinckt 
— [/cosat + — cosaj) + BsinatJsQ. (7) 

where the rotation angle is at — a{t) — \ \d\\{t — to) and 
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B = 



2 sin^ 9 cos^ (f> sin^ 9 sin(20) sin(26') cos ( 
sin^ 6* sin(20) 2 sin^ sin^ (j!) sin(20)sin(; 
sin(20) cos (J) sin(20) sin (j) 2 cos^ 9 

cos 9 — sin 9 sin c 

— cos 9 sin 9 cos 

sin 9 sin 6 — sin 9 cos (h 



and the unit vector d is expressed in spherical coordi- 
nates, d = (sin^cos0, sin^sin^, cos(^)^. 

Without loss of generality we shall assume that we can 
initialize the system in the state Sq = |0) = (0, 0, 1)^ with 
respect to the measurement basis, and that we can exper- 
imentally determine the value of the z-component of the 




FIG. 2: Bloch sphere and arrangement of rotation axes. We 
can define the relative angular displacement of the vectors dk 
representing applied Hamiltonians with respect to a reference 
vector dr (T^z-axis). 



Bloch vector z — (az) by repeated measurements. Insert- 
ing the expressions for A and B into Eq.(jZj), a rotation 
about the axis d by the angle at = \ \d\\{t — to) transforms 
the initial state Sq into s(at) whose z-component is 



z(at) — cos{at) sin^ 9 + 



cos 



(8) 



= 0. 

=.2 , 



Otherwise, it oscillates 
2 , 



z[at) is constant exactly if 

and with minimum Zmin — cos^ 9 — sin'^ 9 — cos(20) for 
cos(Q;t) — —1 or at = (2n -|- l)7r (for some integer n). 
Hence, ii 9 ^ then we can experimentally find ||d|| and 
9 by determining the tmin = pyj and (Fig. P). 

For a single rotation it is sufficient to determine ||d|| 
and 9 and set cj) = 0. For multiple rotations about differ- 
ent axes dfe, however, we must also determine the relative 
azimuthal angles 4>k = 4>'k~ 't'r with respect to a fixed ref- 
erence axis dr (Fig. |2) . The reference axis should have a 

^] and we shall focus on this 
ving determined ||d|| and 9 of 
the axes, and chosen a suitable reference axis d^, we can 
initialize the system in the state Si — {cos/3, sin/?, 0)"^ by 

rotating Sp about the axis d^ by — cos~^ ( cos(2e''j-i 

where P = tan^-'^ (— -^Z— 2 cos(26'r) sec 9r)- Rotating Si by 
various angles a about the axis d^ then gives z{a) — 
C(l — cos a) -I- _Dsina with C = ^ sin(26'fc) cos((/)fc — j3) 
and D = sin^^ sin(0fc — /?), or equivalently 



vertical tilt angle 9r_^ 
case in this paper [ij 



z{a) = 7 sin(a -f- (5) — 7 sin 5, 



(9) 



where C — 7 sin 5 and D = 7COS(5. Hence, we can obtain 
(j) from experimental data for z{a) by finding 5 and 7. 

To determine the evolution of the system subject to the 
several external fields we choose several field strengths 
fm = 1, ■ ■ ■ ,L) for each control /,„ (m = 1, . . . , M) 
and compute the rotation axes 



by finding 9^^]^^ and ||dg'^^|| for all control settings, 
choosing a reference axis, and determining all the rel- 



do 



(10) 



ative azimuthal angles 



-'O+m 



using the strategy outlined 



above. Then we plot the x, y and z-components of d^^^ 

for £ = 1,2, ... ,L versus the field strengths fm for each 
field, and fit a straight line to each set of data points. 
The vertical axis intercepts of the lines then give the x, 
y, and z-component of do, and their slopes determine the 
X, y, and z-component of d,„. 

To evaluate our strategy and find the best ways to ex- 
tract the information from noisy data we use computer 
simulations. We choose various sets of Hamiltonians, use 
the proposed strategy to identify them and compare these 
results with the actual values. Experiments are simu- 
lated by computing individual measurements of s{t) 
according to Eq.(|7|), generating N pseudo-random num- 
bers r„ g [Oil]j where N is the number of times the 
experiment is repeated, and taking the result of the n}^ 
measurement M„ to be 1 if r„ < (1 -I- z)/2 and other- 
wise. To account for measurement errors we introduce a 
symmetric error probability rj 6 [0, 1] indicating the fre- 
quency of erroneous measurement results, i.e., register- 
ing 1 when we should have measured 0, and vice versa. 
This is achieved by choosing another N random numbers 
e„ G [0, 1], and changing Af„ to 1 — M„ whenever e„ < rj. 
Finally, we set (cr^> = -1 + ;| Y.n= 



We illustrate our Hamiltonian identification strategy 
in detail with a specific test system: 



2Ho 
2H2 



0.2cr^ 

0.2cr^ 



0.1a, 
0.9a, 



2Hi = a^ + Q.%ay + 0.la, 



(11) 

This Hamiltonian may arise, for instance, for a charge 
qubit j]j consisting of two quantum dots sharing a single 
electron, and two voltage gates B and S, intended to en- 
able us to change the potential barrier between them and 
induce asymmetries in the double well potential. The 
measurement basis states |0) and |1) are the localized 
ground states of the charge in either of the two dots, and 
projective measurements onto them can be performed via 
a single electron transistor placed next to one of the dots, 
-ffo indicates some tunneling between the dots when no 
control voltages are applied to the gates and a slight 
potential asymmetry, which might be due to imperfect 
placement of the donor impurities. Hi induces tunnel- 
ing between the dots as desired but also induces a slight 
potential asymmetry, which might occur if the electrode 
was placed slightly off-center. H2 produces the desired 
energy level shift but also increases tunneling slightly. 

The first step towards identifying the Hamiltonians Hq, 
Hi and H2 involves finding the rotation frequencies ||do|| 



as well as lid, 



0+1 1 



^0+21 



and the angles 9^, OqL, 9, 



for several gate voltages and . We do this by 
applying each gate voltage for various periods of time 
tj and finding z{tj) = {<Jz{tj)) (Fig- EI- An estimate of 
the rotation frequency lo is obtained by taking the dis- 
crete Fourier transform (DFT) of z{tj) — ■jJ2j=i ^i^j) 
and identifying its maximum. We then estimate the 
time = ^ when z{at) is expected to assume its first 
minimum, acquire (circa 150) additional data points for 
t £ [0.85t,, 1.15^*] and determine the minimum Zmin by 
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FIG. 3: Identification of ||d|| and 61 for d = do + /j'^'M 



Shown is simulated data for d. 



(10) 

0+2 



for /^"'^ = 0.5. The top 



graph is used to estimate the rotation frequency using the 
DFT, the bottom left plot shows the Fourier transform, and 
the bottom right plot the magnification of the region of inter- 
est as well as the parabola that provides the best least-squares 
fit to the data. 




Angle a of 2nd rotation about d*"?] 



FIG. 4: Identification of for d = do -I- /i^°''d2. The points 
*, are a coarse sampling of the curve from which the first 
axis crossing can be estimated (red circle). This then gives 
estimates of the turning points and a^i^^ which can be 

refined by re-sampling (dots) in their vicinity and fitting a 
parabola (red curves). From this we determine the vertices 



i) and (^Oimax,^max 



fitting a parabola to these data points; 9 = ^ cos^^{zjnin)- 
This proved to be considerably more accurate for noisy 
data than fitting a function of the form Eq.Q to the 
data points and reading off the frequency and minimum 
from this curve directly. 



Next we determine the angles (j>Q^j^ 
reference axis, chosen to be dg here. 
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In all of the 
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FIG. 5; Identification of the components of do and dm- 
Typical simulated experimental run for the system defined 
by Ea. pi^ and the estimated components of d for different 
field settings. Shown are the best straight line fit (dashed 
lines) and the actual values (solid line). Note the excel- 
lent agreement between the actual values and the best fit 
obtained. A typical run of the program yielded: do — 
(0.1987, 0.0064, 0.0992)^ di = (0.9859, 0.9122, 0.1149)^, da = 
(0.2081,0.0170,0.8957)^, i.e., the Hilbert-Schmidt norm er- 
rors ||d^*-d^^*l| were less than 3% (0.0066, 0.0238 and 0.0193 
for m = 0, 1, 2 respectively). 



following experiments the system is initialized in state 
Si = (cos/?,sin/3,0)-^ by applying a suitable rotation 



about the reference axis. In our case the state Sq is 
rotated to Si by letting it evolve freely for « 8.1379 
time units yielding f3 « —1.0489. We then apply each 
field /m"* for various times to achieve rotations by vari- 
ous angles a (Fig. 0]). This yields the parameters 7 — 
Zm.ir,.) and (5 = TT - \{amin + ctmax) in Eq.lP) 
e.g., by setting D = 7 cos (5 



2 (^maa: ^min) and 5 

from which we can obtain 



and (j) = —(3 — svoT (£'/sin(6')) where is the vertical 
tilt angle of the rotation axis determined in the previous 
step. 

Finally, having determined all the relevant parameters 
of the rotation axes, we convert the data to Cartesian co- 
ordinates, and plot the x, y and z components of dp^^^ 

for each field as a function of /m"* (Fig. EJ . From the 
straight line fitting, we can then estimate the Hamiltoni- 
ans. 

In conclusion, we have presented a stepwise procedure 
to determine the dynamical response of a Hamiltonian 
2-level quantum system to control fields when we only 
have access to measurement in a single fixed basis. This 
overcomes a weakness in the assumptions of QPT and 
thus enables us to bootstrap the procedure to estimate 
parameters for systems consisting of arrays of qubits. Nu- 
merical results indicate that the procedure is fairly robust 
for Hamiltonians of physical interest. Extensions of the 
procedure higher-dimensional or dissipative systems may 
be possible. 
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